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^ , The confinement problem has been solved in the anisotropic (2+ l)-dimensional SU(A^) 

^ ! Yang-Mills theory at weak coupling. In this paper, we find the low-lying spectrum for 

I N = 2. The lightest excitations are pairs of fundamental particles of the (1 + 1)- 

O ■ dimensional SU(2) x SU(2) principal chiral nonlinear sigma model bound in a linear 
potential, with a specified matching condition where the particles overlap. This match- 

^ . ing condition can be determined from the exactly-known S-matrix for the sigma model. 
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1 Introduction 



In recent papers, some new techniques have been developed for calculating quantities in 
a (2+l)-dimensional SU(A^) gauge theories [1], [2], [3]. These techniques exploit the fact 
that in an anisotropic limit of small coupling, the gauge theory becomes a collection of 
completely-integrable quantum field theories, namely SU(A^) x SU(iV) principal chiral 
nonlinear sigma models. These integrable systems are decoupled, save for a constraint 
which is necessary for complete gauge invariance. In the case of = 2, is possible to 
perturb away from integrability, using exactly-known off-shell matrix elements of the 
integrable theory. 

Though the gauge theory we consider s not spatially-rotation invariant, it has fea- 
tures one expects of real (3-|-l)-dimensional QCD; it is asymptotically free and confines 
quarks at weak coupling. Thus the limit of no regularization is accessible. 

One can formally remove the regulator in strong-coupling expansions of (2 + 1)- 
dimensional gauge theories; the vacuum state in this expansion yields a string tension 
and a mass gap which have formal continuum limits. This can be done in a Hamiltonian 
lattice formalism or with an ingenious choice of degrees of freedom and point- 
splitting regularization [5] . This leaves open the question of whether these expressions 
can be trusted at weak coupling (more discussion of this issue can be found in the 
introduction of reference [2]). In particular, one would like to rule out a deconfinement 
transition, or different dependence of physical quantities on the coupling (as in compact 
QED There is a proposal for the vacuum state [7], in the formulation of reference 
[5] which seems to give correct values for some glueball masses [B|, but this proposal 
evidently requires more mathematical justification. 

In this paper, we will work out the masses of the lightest glueballs for the case 
of gauge group SU(2). Our method would also work in principle for SU(A^) gauge 
theories, and our reason for choosing iV = 2 is that the analysis is simplest for that 
case. 

The basic connection between the gauge theory and integrable systems is most easily 
seen in axial gauge [1] . The string tensions in the x^-direction and x^-direction (which 
we called the horizontal and vertical string tensions, respectively) for very small (^q, 
were found by simple physical arguments. The result for the horizontal string tension 
was confirmed for gauge group SU(2), and additional corrections in q'q were found [2], 
through the use of exact form factors for the currents of the sigma model. String 
tensions for higher representations can also be worked out, and adjoint sources are not 
confined [3]. 

Careful derivations of the connection between the gauge theory and integrable sys- 
tems use the Kogut-Susskind lattice formalism [I], [2]- A shorter derivation was given 
in reference [5], which we summarize again here. The formalism is essentially that of 
"deconstruction" fiO\. 

The Yang-Mills action is / d^C, where the Lagrangian is £ = ^^TyFq^ + ^TtFq2 — 
2^TrFf2, and where Ao,Ai and Ai are SU(A^)-Lie-algebra- valued components of the 
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gauge field, and the field strength is F^^, = d^A^, — d^A^ — i[A^, Ai,]. This action 
is invariant under the gauge transformation Afj_{x) —>■ ig{x)^^[d^ — iAf^{x)]g{x), where 
g{x) is an SU(A^)-valued scalar field. We take e' 7^ e, thereby losing rotation invariance. 

We discretize the 2-direction, so that the takes on the values = a, 2a, 3a ... , 
where a is a lattice spacing. All fields are considered functions of We 
define the unit vector 2 = (0,0, 1). We replace ^2(0;) by a field U{x) lying in SU(A^), 
via U{x) ~ exp — my42(a;). There is a natural discrete covariant-derivative operator: 
V^^{x) = d^ii{x) - iA^{x)ii{x) + iil{x)A^{x + 2a), /i = 0, 1, for any N x N complex 
matrix field il(x). The action is 5 = / dx'^ J dx^ J2x^ ^ ^ where 

^ = :v^TrFo\ + ^T^\V^U{x)\^V^U{x) - ^T^\V^U{x)\^V^U{x) , (1.1) 

and where g\ = e^a and ((^q)^ = ^' ^'j^- The Lagrangian (11.11) is invariant under the gauge 
transformation: A^{x) ig{x)~^[dfj, — iAfj_{x)]g{x) and U{x) — > g{x)~^U{x)g{x + 2a) 
where again, g{x) G SU(iV) and /i is restricted to or 1. The bare coupling constants go 
and g'f) are dimensionless. We recover from (11.11) the anisotropic continuum action in the 
limit a — > 0. The sigma model field is U{x^, x^,x'^), and each discrete x'^ corresponds to 
a different sigma model. The system (11.11) is a collection of parallel (1 + l)-dimensional 
SU(A^) X SU(A^) sigma models, each of which couples to the auxiliary fields Aq, Ai. 
The sigma- model self- interaction is the dimensionless number g^. 

We feel it worth commenting on the nature of the anisotropic regime and how it 
is different from the standard (2 -|- l)-dimensional Yang-Mills theory. The point where 
the regulator can be removed in the theory is gQ = go = 0. This point can be reached 
in our treatment, but only if 

[g'^y < le-M(^;gA^) . (1.2) 
fl'o 

The left-hand side and ride-hand side are proportional to the two energy scales in the 
theory (the latter comes from the two- loop beta function of the sigma model). Thus 
our method cannot accommodate fixing the ratio g'o/go, which is natural in standard 
perturbation theory ^Jj. This is why the mass gap is not of order e, e' and the string 
tension is not of order e^, (e')^. 

We now discuss the Hamiltonian in the axial gauge Ai = 0. The left-handed and 
right-handed currents are, j^(a;)b = iTTthd^U{x) U{xy a.ndj^{x)b = iTTthU{xydfj_U{x), 
respectively, where /i = 0, 1. The Hamiltonian obtained from (II. ip is Hq + Hi, where 

Ho = J2 I dx'^,{[j^o{^\f + \j\{x\f} , (1.3) 
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where — $fc = t is the temporal gauge field, and where in the last term we have 
inserted two color charges - a quark with charge q at site u and an anti-quark with 
charge q' at site v. Some gauge invariance remains after the axial- gauge fixing, namely 
that for each 



where Q{x'^)b is the total color charge from quarks at and \E' is any physical state. To 
derive the constraint (11. 5p more precisely, we started with open boundary conditions 
in the 1-direction and periodic boundary conditions in the 2-direction, meaning that 
the two-dimensional space is a cylinder |1], [2]- 

From (11.41) we see that the left-handed charge of the sigma model at x^ is coupled 
to the electrostatic potential $, at x"^. The right-handed charge of the sigma model 
is coupled to the electrostatic potential at + a. The excitations of Hq, which we 
call Fadeev-Zamoldochikov or FZ particles, behave like solitons, though they do not 
correspond to classical configurations. Some of these FZ particles are elementary and 
others are bound states of the elementary FZ particles. An elementary FZ particle has 
an adjoint charge and mass mi. An elementary one-FZ-particle state is a superposition 
of color-dipole states, with a quark (anti-quark) charge at x^,x^ and an anti-quark 
(quark) charge at x^,x^ + a. The interaction Hi produces a linear potential between 
color charges with the same value of x^. Residual gauge invariance (11.51) requires that 
at each value of x^, the total color charge is zero. If there are no quarks with coordinate 
x^, the total right-handed charge of FZ particles in the sigma model at x^ — a is equal 
to the total left-handed charge of FZ particles in the sigma model at x^. 

The particles of the principal chiral sigma model carry a quantum number r, with 
the values r = 1, . . . , — 1 [21]. Each particle of label r has an antiparticle of the 
same mass, with label N — r. The masses are given by 

rrir = mi— — ^, mi = KA(gQN)~^^'^e + non — universal corrections , (1.6) 



where A' is a non-universal constant and A is the ultraviolet cut-off of the sigma model. 

Lorentz invariance in each x°,x^ plane is manifest. For this reason, the linear 
potential is not the only effect of Hi. The interaction creates and destroys pairs of 




(1.5) 




3 



elementary FZ particles. This effect is quite small, provided that q'q is small enough. 
Specifically, this means that the square of the 1 + 1 string tension in the x^-direction 
coming from Hi is small compared to the square of the mass of fundamental FZ particle; 
this is just the condition (11.21) . The effect is important, however, in that it is responsible 
for the correction to the horizontal string discussed in the next paragraph in equation 

dLHD. 

Simple arguments readily show that at leading order in g^^, the vertical and hori- 
zontal string tensions are given by 

a la^ 

respectively, where Cat is the smallest eigenvalue of the Casimir of SU(A^). These naive 
results for the string tension have further corrections in g'^^ which were determined for 
the horizontal string tension for SU(2) [2]: 



<T„ = - I - 



-1 

(1.^ 
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The leading term agrees with (11.71) . This calculation was done using the exact form 
factor for sigma model currents obtained by Karowski and Weisz [12] . The form factor 
can also be used to find corrections of order ((^g)^ to the vertical string tension; this 
problem should be solved soon. If the reader is not familiar with form-factor techniques 
in relativistic integrable field theories, a self-contained review is in the appendix of 
reference 0. 

Another recent application of exact form factors to the (2 + l)-dimensional SU(2) 
gauge theory is reference pJJ , in which form factors of the two-dimensional Ising model 
[13] are used to find the profile of the electric string near the high-temperature decon- 
fining transition, assuming the Svetitsky-Yaffe conjecture [T5]. 

A rough picture of a gauge-invariant state for the gauge group SU(2) with no quarks 
is given in Figure 1. For N > 2, there are more complicated ways in which strings can 
join particles. For example, a junction of strings is possible. Figure 1 is inaccurate 
in an important respect; the "ring" of particles held together by horizontal strings is 
extremely broad in extent in the x^-direction compared to the x^-direction. This is 
because cth -C cry- 

The hghtest states have the smallest number of particles, by virtue of cth <^ cry- 
Thus the lightest glueballs are pairs of FZ particles with the same value of x^. For 
small enough g'^, the very lightest state has a mass well- approximated by 2mi. The 
purpose of this paper is to find the leading corrections in {g'^Y this result. This 
will be done using the S-matrix of the sigma model and the WKB formula. There 
are further small corrections, due to the softening of the potential near where particles 
overlap, which we do do not determine. 

It is clear that the lightest bound states of FZ particles are (1 + l)-dimensional in 
character. If we formulated a gauge theory in which was fixed in f/(a;°, x^, x^), we 
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would find the same spectrum, as a function of mi and cth- In the Kogut-Susskind 
lattice formulation, a long row of plaquettes with open boundary conditions is a regular- 
ized gauge theory of this type. The only real difference between this (1 + 1) dimensional 
model and that we study is that uh will receive different corrections of order ((^q)^- 




Figure 1. A glueball state is a collection of heavy particles, held weakly together 
by strings. The horizontal coordinate is and the vertical coordinate is x"^. 

In the next section we will discuss the wave function of an unbound pair of FZ 
particles. We find that this is described by phase shift for the color-singlet sector. In 
Section 3, we determine the bound-state spectrum. The problem we solve is very similar 
to that of two particle-states of the two-dimensional Ising model with an external 
magnetic field \i6\ (for a good summary of this problem, see reference |IT]); the only 
genuine difference is the presence of a matching condition where the particles overlap. 
This matching condition comes from the phase shift of the scattering problem. We 
present our conclusions in Section 4. 

2 Scattering states of FZ particles 

The lightest glueball state, as discussed above, is simply a pair of FZ particles located 
at the points (x^,x^) and {y^,x^) and bound in a linear potential. Residual gauge 
invariance (11. 5p . demands that the state be a color singlet. To begin with, however, 
we simply write the form of a free state of two particles. 

The state of the SU(2) x SU(2) ~ 0(4) nonlinear sigma model with a particles of 
momenta pi and p2 and quantum numbers ji and j2 (which take the values 1,2,3,4) 
is described by the wave function 



'^piP2 {-^ J y ) 




x^ < 
x^ > 



(2.1) 
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where Aj-^j^ is an arbitrary set of complex numbers and Sj^j^{pi,p2) is the two-particle 
S-matrix. We have not yet imposed (11 .SI) . 

The wave function (12.11) is written in a form where the 0(4) symmetry is manifest. 
It is straightforward to write it in a form where the left SU(2)l and the right SU(2)pi, 
symmetries are manifest, by writing 

ij,,,,{x\ y^S = E -^(5^^.^ - ^O^iSlf - iaiir ^,,,,{x\ y%,, (2.2) 

describing a pair of color dipoles, one with quantum numbers a, b and the other with 
quantum numbers c, d, where , j = 1, 2, 3 denotes the Pauli matrices. 

We impose the physical state condition (II. 5p on (12. 2p by requiring that a = b and 
c = d and summing over these colors. The projected wave function is, up to an overall 
constant, 

^p,p,{x\y') = I ^ ^i^yi , (2.3) 

where 5*0(^1,^2) is the singlet projection of the 0(4) S-matrix. This S-matrix was first 
obtained by Zamolodchikov and Zamolodchikov [18J . A useful form is given in reference 

So{PuP2) = So{e) = — expi / — — — sm— , (2.4) 

TT + It^ Jq 4 1 + e« TT 

where the relative rapidity 6 is given hj 6 = 62 — 61, pi = m sinh 61, p2 = m sinh 62 and 
where we denote the particle mass mi, given by (11.61) . by m (because there is only one 
mass for the case of = 2). A derivation of (12.41) is in the appendix of reference [2]. 

The singlet S-matrix is just given by a phase shift (p{9): Sq{9) = expi(f){9). The 
phase shift has a simple form in the low-energy, non-relativistic limit, \pi — P2I ttl. 
In this limit, 6 ^ \pi — p2\/m. The integral on the right-hand side of (12.41) can be done 
by Taylor expanding in \pi —p2\/m yielding 

3-21n2, , ^/b-r|2\ ^ ^ 

= 0(pi,P2) = 7r \pi-P2\ +0r-^^-^) . 2.5 

Txm V Tn"^ J 



3 The low-lying glueball spectrum 

Let us now consider the states of a bound pair of FZ particles in the potential V{x^, y^) = 
2a}i\x^ — y^\ (the reason for the factor of two is simply that the particles are joined by 
a pair of strings). We use the non-relativistic approximation, used to find (12.51) . For 
our problem, the horizontal string tension times the size of a typical bound state is 
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small compared to the mass, by fll.2p . This justifies the non-relativistic approximation 
for low-lying states. The mass of a low-lying glueball is given by 

M = 2m + E , 

where E is the energy eigenstate of the two-particle problem. 

Let us introduce center-of-mass coordinates, X = {x^ + y^)/2 and x = — x^. The 
reduced mass of the system is m/2. We factor out the phase depending on X, leaving 
us only with a wave function depending on x. The Schrodinger equation we consider 
is 

m dx'^ 

with a matching condition at x = between the wave function 'il){x) at x > and the 
wave function at x < 0. There is actually a further complication, which we do not 
consider here; the potential changes slightly in the region where x ~ 0. This is due 
to the fact that the color charge is slightly smeared out. This smearing out can be 
calculated from the form factor p2] . 

Our results (12. 3p . (12. 5p for the unbound two-particle state, tell us that for x^ ~ y^, 
where the effect of the potential can be ignored, the bound-state wave function in the 
center-of-mass frame will be of the form 

, ■ ( cos{px + uj) , X < . , 

= { r , A/ \1 n ; (3-1) 

COS[— px + 00 — (p[p)\ , X > U 

for some angle co, where p = pi — P2 and (f){p) = tt — |p| + 0(|pp/m^). The value 

of p near x = is given hj p = (mEY^'^, where E is the energy eigenvalue of the state. 
This is the matching condition between the wave function for x > and for x < 0. 

The wave function for x < an Airy function. So is the wave function for x > 0. 
We therefore obtain the approximate WKB form 



■ip{x) 



3(2maH)'/'(x^ 2^ 
-f(2maH)i/^(^-x)3/2 + | 



X < 

, (3.2) 

X > 



for some constants C and C. The expression (13.21) can be made to agree with (13.11) for 
small X, provided the generalization of the Bohr-Sommerfeld quantization condition 

2(m)i/2 . 3-21n2 f l\ , ^ 

3(Th Trm^/^ \ 2/ 

is satisfied hj E = En- The only new feature in this semi-classical formula is the 
second term, produced by the phase shift. Absorbing the horizontal string tension in 
the energy, by defining m„ = Encr^^^^ , this cubic equation becomes 

2(m)^ 3-21n2 V3^V2 _ + 1^ , - Q 
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The second term can be ignored for sufficiently small o"h, i.e. sufficiently small Qq. 

There is a unique real solution of the cubic equation fl3.3p for a given integer n > 0, 
because 3 — 21n2 = 1.613706 > 0. The low-lying glueball masses are given by 

a/3 _ 3(3-21n2)ffH 1/3I ^ 



M„ = 2m + En = 2m + 



where 



4^1/2 



+ 



n 



1 

+ 7^ 



3(3-21n2)aH 



2Tcm 



1/2 



(3.4) 



(3.5) 



4 Conclusions 

We have identified the low-lying glueballs of the anisotropic Yang-Mills theory in (2 + 
1) dimensions as bound pairs of the fundamental massive particles of the principal 
chiral nonlinear sigma model. We found a matching condition for the bound-state 
wave function at the origin, which when combined with elementary methods yields the 
spectrum of the lightest states. 

There are other aspects of the two-particle bound-state problem we have not con- 
sidered here. First, the potential is not precisely linear in the region where the two 
particles are close together. The corrections to the potential can be determined us- 
ing form factors. This will slightly modify (13. 3p . A completely different issue is that 
there are small corrections to the form factors themselves, coming from the presence 
of bound states. This, in turn, will give a further correction to the horizontal string 
tension found in [2]. Such corrections to form factors in theories close to integrability 
were ffist discussed by Delfino, Mussardo and Simonetti [12]. The bound-state energies 
proliferate between 2m and 4m, as q'q — > 0. Our method breaks down as the bound- 
state mass reaches 4m, because the bound state develops an instability towards fission 
into a pair of two-particle bound states. This is analogous to the situation for the Ising 
model in a field [TB], [I?] as we stated earlier. It should be worthwhile to understand 
the relativistic corrections to the bound-state formula, along the lines of the work of 
Fonseca and Zamolodchikov [20]. 

A similar calculation is possible for SU(A^). The exact S-matrix of the principal 
chiral nonlinear sigma model is known for > 2 [21] . An interesting feature is that the 
phase shift should vanish as A^ — > 00, with g^N fixed, meaning that the wave function 
would be continuous where FZ particles overlap. 

It would be interesting to study the scattering of a glueball by an external particle. 
If the scattering is sufficiently short range, the FZ particles could be liberated from the 
glueball, after which hadronization would ensue. 

The results of this paper and of references [I] and [^ may be extendable to the 
standard (2 + l)-dimensional isotropic Yang- Mills theory with Qq = qq. The strat- 
egy we have in mind is an anisotropic renormalization procedure. At the start is a 
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standard field theory with an isotropic cut-off. By anisotropically integrating out high- 
momentum degrees of freedom, the isotropic theory will flow to an anisotropic theory 
with a small momentum cut-off in the x^-direction and a large momentum cut-off in the 
direction. If the renormalized couplings satisfy the condition (II. 2p . we could apply 
our techniques. A check of such a method would be approximate rotational invariance 
of the string tension. This would give an analytic first-principles method of solving 
the isotropic gauge theory with fixed dimensionful coupling constant e, and no cut-off. 
The only other analytic weak-coupling argument for a mass gap and confinement in 
(2 + l)-dimensions, namely that of orbit-space distance estimates, discussed by Feyn- 
man [22], by Karabali and Nair in the second of references f5], and by Semenoff and 
the author [23] is suggestive, but has not yielded definite results yelQ. 
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